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Abstract 

An  integrated  procedure  for  math  modeling  and  power  control  strategy  design  for  a  fuel  cell  hybrid  vehicle  (FCHV)  is  presented  in  this  paper. 
Dynamic  math  model  of  the  powertrain  is  constructed  firstly,  which  includes  four  modules:  fuel  cell  engine,  DC/DC  inverter,  motor-driver,  and 
power  battery.  Based  on  the  mathematic  model,  a  power  control  principle  is  designed,  which  uses  full- states  closed-loop  feedback  algorithm.  To 
implement  full-states  feedback,  a  Luenberger  state  observer  is  designed  to  estimate  open  circuit  voltage  (OCV)  of  the  battery,  which  make  the 
control  principle  not  sensitive  to  the  battery  SOC  (state  of  charge)  estimated  error.  Full- states  feedback  controller  is  then  designed  through  analyzing 
step  responding  of  the  powertrain  and  test  data.  At  last  of  the  paper,  the  results  of  simulation  and  field  test  are  illustrated.  The  results  show  that 
the  power  control  strategy  designed  takes  into  account  the  performance  and  economy  characteristics  of  components  of  the  FCHV  powertrain  and 
achieves  the  control  object  excellently. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

In  recent  10  years,  with  the  problems  of  energy  shortage  and 
air-pollution  become  more  and  more  serious,  interest  in  alterna¬ 
tive  automotive  powertrain  increases  steadily.  As  one  substitute 
for  traditional  vehicle,  fuel  cell  vehicles  (FCV)  have  become  a 
research  hotspot  in  area  of  clean  vehicle  technology.  The  main 
purpose  of  this  paper  is  to  describe  the  model-based  development 
processes  of  the  FCV  power  control  strategy  which  is  applied 
in  Tongji  University  START3  prototype  fuel  cell  hybrid  vehicle 
(FCHV). 

As  the  automotive  engine,  fuel  cells  have  some  problems  such 
as  slow  responding,  difficulty  cold  starting,  and  braking  energy 
cannot  be  recycled.  So  power  battery  is  applied  in  FCV  pow¬ 
ertrain  usually.  This  is  so-called  fuel  cell  hybrid  vehicle.  Most 
major  vehicle  manufactures  have  developed  prototype  FCV  for 
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technology  evaluation  and  demonstration  purpose  [1-3].  How¬ 
ever,  since  these  vehicles  are  not  meant  for  production,  they 
are  put  together  swiftly.  Model-based,  systematic  development 
processes  that  are  commonly  adopted  for  traditional  (internal 
combustion  engine  powered)  vehicles,  are  rarely  explored  to 
their  full  potential  in  the  development  of  these  prototype  vehi¬ 
cles. 

Academic  research  groups  are  also  working  actively  on  this 
emerging  technology.  Guezennec  et  al.  [4]  solved  the  supervi¬ 
sory  control  problem  of  a  FCHV  as  a  quasi- static  optimization 
problem  and  found  that  hybridization  can  significantly  improve 
the  fuel  economy  of  FCVs.  Mohsen  Mohammadian  et  al.  [5]  deal 
a  control  strategy  developed  for  optimizing  the  energy  flow  by 
using  evolutionary  algorithms  implemented  on  a  fuel  cell  hybrid 
vehicle  to  reach  the  best  performance,  fuel  economy,  emission 
and  acceptable  operation  of  this  hybrid  structure.  The  National 
Renewable  Energy  Laboratory  included  a  fuel  cell  system  option 
in  its  simulation  model  ADVISOR  [6] ,  which  can  simulate  FCV s 
and  HE  Vs  with  ease. 

Up  to  today,  few  papers  have  published  detailed  models  of 
a  FCHV  powertrain,  and  model-based,  systematic  development 
processes  was  not  used  in  control  design  for  FCHV  powertrain. 


822 


Y.  Haitao  et  al.  /  Journal  of  Power  Sources  180  (2008)  821-829 


Fig.  1.  Schematic  configuration  of  the  START3  FCHV. 


The  major  feature  of  this  study  is  model-based  control  system 
design  for  FCHV  powertrain.  Which  is  carried  out  by  three  steps: 
firstly,  a  state-space  mathematic  model  of  the  FCHV  powertrain 
is  formulated  based  on  its  equivalent  circuit  model  and  test  data. 
Then,  power  control  strategy  is  developed  based  on  the  mathe¬ 
matic  model,  which  uses  full-states  feedback  algorithm.  At  last, 
the  strategy  is  validated  by  simulation  and  field  test. 

2.  State-space  mathematic  model  of  the  powertrain 

2.7.  Equivalent  circuit  model  of  the  powertrain 

The  FCHV  studied  in  this  paper  is  Tongji  University  START3 
prototype  vehicle.  The  powertrain  configuration  of  START3  is 
shown  in  Fig.  1.  A  permanent  magnet  synchronization  motor 
drives  wheels  through  a  slow-down  gear-box.  The  electric  motor 
inverter  uses  vector  control  strategy  and  the  electric  motor  can 
work  in  four  quadrants.  A  DC-DC  converter  (DC/DC)  is  placed 
between  fuel  cell  engine  (FCE)  and  the  motor  inverter,  its  pur¬ 
pose  is  to  compensate  the  over  soft  output  characteristic  of  FCE, 
namely  resistance  matching.  A  Li-ion  battery  is  parallel  con¬ 
nected  to  main  circuit,  as  assistant  power  source,  its  function 
is  recycling  braking  energy,  providing  peak  power,  and  helping 
with  the  transient  part  of  power  load. 

In  the  powertrain  of  START3,  the  DC/DC  is  current  control 
and  the  motor  is  torque  control,  so  both  of  them  can  be  regarded 
as  a  constant  current  source.  Bus  voltage  of  the  main  circuit 
is  equivalent  to  output  voltage  of  battery.  An  equivalent  circuit 
model  can  be  abstracted  consequently,  as  shown  in  Fig.  2.  In 
this  equivalent  circuit  model,  battery  is  represented  by  a  third- 
order  equivalent  circuit  model  which  is  proposed  and  validated 
in  [7]. 

2.2.  Testbench 

Most  components  of  START3  powertrain  were  designed 
specifically.  They  were  bench-tested  for  factory  acceptance 
before  vehicle  integration.  A  test  bench  was  constructed,  which 
enables  more  precise  and  controlled  testing  of  individual 
components.  In  addition  to  subsystem  function  verifications, 
special-purpose  tests  were  also  designed  to  obtain  the  neces¬ 
sary  data  to  build  dynamic  model  and  efficiency  map  of  the 
subsystems. 


The  component  configuration  of  test  bench  for  fuel  cell 
engine  and  electric  motor  is  roughly  same  as  the  powertrain 
of  START3  while  an  electric  power  load  system  (supplied  by 
Arbin)  was  used  to  simulate  the  power  load  or  power  source. 
The  DC/DC  and  battery  pack  were  tested  using  the  test  bench 
which  was  designed  specially. 

2.3.  Mathematic  modeling 

2.3.1.  FCE  model 

In  the  powertrain  control  system,  the  input  command  of  FCE 
is  power  setting  Ffseb  the  output  is  upper  limit  of  output  power 
T^Lmt-  One  fact  should  be  known  is  that  in  the  powertrain  sys¬ 
tem,  T^Lmt  only  denote  the  real-time  output  capacity  of  FCE, 
but  not  the  real-time  output  power  of  FCE.  Real-time  output 
power  of  FCE  is  influenced  by  power  load  and  battery  SOC 
(state  of  charge).  According  to  [8],  FCE  is  a  nonlinear  dynamic 
system,  and  in  fuel  cell  control  level  FCE  model  must  be  non¬ 
linear  and  very  complicated.  But  in  powertrain  control  level  too 
complicated  FCE  model  is  unnecessary.  In  the  powertrain  con¬ 
trol  system,  the  relationship  between  FTmt  and  Pfset  is  just  the 
FCE  model.  According  to  [8]  and  bench  test  data  (Fig.  3),  the 
relationship  between  them  can  be  simplified  as  a  first-order  sys¬ 
tem  approximately,  as  the  transfer  function  shown  in  Eq.  (1).  In 
Fig.  3,  it  also  can  be  seen  that  time  constant  of  the  first-order 
function  is  just  little  different  in  anterior  time  and  posterior  time. 


Fig.  2.  Equivalent  circuit  model  of  the  FCHV  powertrain.  Where  Im:  current 
of  the  motor;  4cf.‘  current  of  the  DC/DC;  4:  current  of  the  battery;  Ro,Ri,R2i 
equivalent  resistant;  Cq,Ci  El'-  equivalent  capacitance;  E:  equivalent  voltage. 
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Fig.  3.  Relation  between  Pfset  and  /Ymt  in  bench-test  (where  Pfset  =  40  Kw). 


the  FCE  and  DC/DC  are  looked  as  one  current  source.  According 
to  Eq.  (1)  and  (2)  this  current  source  can  be  denoted  by  transfer 
function,  as  shown  in  the  following  equation: 


Gxf(s)  = 


Clef1 (s ) 

ClcfSetC?) 


(  Kdcf  \  /  Kfce  \ 
V  TfcfS  +  1  /  v  7fces  +  1  / 


(3) 


where  f/bus'  bus  voltage.  Because  Figs.  3  and  4  show  that 
responding  speed  of  DC/DC  is  much  fast  than  that  of  FCE  (about 
20:1),  Tdcf  can  be  omitted  in  simplified  mathematic  model.  Eq. 
(3)  can  be  simplified  as 


Gjk(s)  = 


Clcf(s) 

ClcfSet(s) 


^dcf^fce 
TfceS  +  1 


T/s+1 


(4) 


where  Kf  =  KdcfKfc&,  Tf=  7fce 


The  simplification  is  advisable  for  mathematic  modeling  and 
control  system  design,  and  it  will  also  be  justified  by  model 
parameters  estimation  and  validation  later. 


^Kfce(s) 


ELmt(s) 

^fSet(s) 


Afce 

7fceS  +  1 


(1) 


where  Kfcc,  7fce:  constants. 


2.3.2.  DC/DC  model 

The  input  command  of  DC/DC  is  current  setting  Idcfseu  the 
output  is  actual  current  Clef.  According  to  the  bench  test  data 
(Fig.  4),  the  relationship  between  them  can  be  simplified  as  a 
first-order  system,  as  the  transfer  function  shown  in  (2).  The 
simplification  is  advisable  for  mathematic  modeling  and  control 
system  design,  and  it  will  also  be  justified  by  model  parameters 
estimation  and  validation  later. 


GKdcf(s)  = 


Clef1 (s ) 

ClcfSet(s) 


^dcf 

TdcfS  +  1 


(2) 


where  Kd cf,  Tdcf:  constants. 

In  practice,  response  speed  of  the  DC/DC  is  slowed  down 
to  protect  FCE,  namely  /dcfset  is  limited  by  /Amt-  To  reflect  this 
relationship  in  mathematic  model,  /dcfset  is  set  to  upper  limit  and 


Fig.  4.  Relation  between  /dcfset  and  /dcf  in  bench-test  (where  /dcfset  =  1 10  A). 


2.3.3.  Motor-driver  model 

The  input  command  of  electric  motor  inverter  is  torque  setting 
Tqset,  the  output  is  actual  torque  Tq/^.  Because  bench- test  data 
shows  that  the  electric  motor  respond  very  fast,  and  the  relation¬ 
ship  between  its  output  and  input  is  linear,  it  can  be  described 
by  a  first-order  system,  as  the  transfer  function  shown  in  the 
following  equation: 


GKmc(s)  = 


TqisOO 

TqsetW 


K 


me 


Tmcs  4“  1 


(5) 


where  Kmc ,  Tmc:  constants.  Because  the  electric  motor  current 
is  function  of  its  torque,  Eq.  (5)  is  modified  as 


G  Km(s)  = 


4n(s) 

Anset(s) 


Tqis(s) 

Tqset(s) 


K 


me 


Tmcs  4”  1 


(6) 


2.3.4.  Battery  model 

As  shown  in  Fig.  2,  the  battery  is  simulated  by  the  third- 
order  RC  circuit  model  which  has  transfer  functions  shown  in 
the  following  equation: 


Gkw  = 


G  Kb\  = 


Gk  b2  = 


Uo(s) 

Ibis) 

Uj(s) 

Us) 

U2(s) 

Us) 


1 

Cos 

Ri 

R\C\s  1 
Ri 

R2C2S  4-  1 


(7) 


where  4:  battery  current,  Co,  U 1,  U2 :  voltages  across  Co,  Ci, 
C2.  According  to  [3],  Uo  is  equivalent  to  open  circuit  voltage 
(OCV)  of  the  battery. 

Uniting  the  Eq.  (4),  (6)  and  (7)  and  taking  inverse  Laplace 
transform,  a  state-pace  model  of  the  powertrain  emerge  as  shown 
in  the  following  equation: 


x  =  Ax  +  Bu 
y  =  Cx  4~  Du 


(8) 


where  state  vector  x  =  [Idcflm  Co  U\  C2]T,  input  vector  u  =  [Efset 
^dcfset  7mse  t]T,  output  vector  y  =  [/dcf  Im  Ubus  Pimt]T.  A,  B,  C,  D 
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are  state  matrix  which  are  defined  in  Eq.  (9). 
1 


- 0  0 

Tf 

0 

1 

— - 0 


0 


1 

Tyne 

1 


B  : 


Co 

1 

Cl 

1 

L  c~2 
o 


Co 

1 

~c[ 

1 

C~2 


0 

0 


1 


RiCi 

0 


— c  0 


Kf 

Tf 

0 


L  me 

0 

0 

0 


Table  1 

Model  parameters  estimation  result 


0 

FAST 

UDDS 

EC_EUDC 

J1015 

mean 

Ee  (%) 

Kfce 

0.9943 

1.0172 

0.9981 

1.0227 

1.008 

9.19 

0 

Tee 

0.9957 

1.0831 

0.9214 

1.1747 

1.0437 

Kdcf 

0.9723 

0.9566 

0.9566 

0.9592 

0.9612 

3.85 

0 

Tdcf 

0.3162 

0.3386 

0.303 

0.3091 

0.3167 

’  Kmc 

0.9697 

0.9852 

0.9814 

0.9863 

0.9806 

3.37 

Tnc 

0.0703 

0.1015 

0.0962 

0.0684 

0.0841 

0 

Co 

896.94 

962.87 

1001.87 

709.88 

892.89 

5.8 

Ro 

0.0993 

0.1058 

0.1067 

0.0909 

0.1007 

1 

Cl 

6.3006 

6.4876 

9.4997 

5.1079 

6.8490 

R2C2 . 

(9)  Ri 

0.0743 

0.074 

0.0601 

0.0385 

0.0617 

c2 

115.63 

127.92 

95.191 

154.26 

123.25 

r2 

0.1771 

0.1683 

0.1777 

0.1134 

0.1591 

C=  [flo  -Ro  1  1  l],D=[0  0] 


2.4.  Model  parameters  estimation 


A  FCHV  powertrain  simulation  model  is  constructed  based 
on  the  mathematic  model  in  Matlab/simulink.  Even  though  each 
subsystem  is  modeled  and  validated  separately,  the  state-space 
model  needs  to  be  checked  to  confirm  that  the  overall  model 
works  properly.  The  validation  is  done  by  comparing  simulation 
result  with  actual  test  data.  The  motor  torque  demand  in  bench- 
test  is  given  to  the  simulation  model  as  the  input. 

To  evaluate  precision  of  the  simulation,  a  variable  is  defined 
in  Eq.  (11),  which  is  similar  to  the  variable  defined  in  Eq.  (10). 


There  are  several  unknown  parameters  in  the  mathematic 
model  developed  above.  To  estimate  these  parameters,  the 
quadratic  least  squares  parameter  estimation  method  is  utilized 
in  this  paper.  It  is  a  method  for  linear  regression  that  determines 
the  values  of  unknown  quantities  in  a  statistical  model  by  mini¬ 
mizing  the  sum  of  the  squared  residuals  (expounded  in  [9]  and 
[10]).  Sample  data  came  from  the  bench- test.  To  evaluate  pre¬ 
cision  of  the  estimation,  a  variable  called  average-error-rate  is 
defined  in  the  following  equation: 


1 

Ee  =  E(e)  =  - 


z(k)  -  hT(k)0 LS 

Z(k)  -  Zbase 

(10) 


where  L:  quantity  of  the  sample  data,  Zbase  is  reference  value  of 
output  in  simulation,  h(k)  is  state  vector  in  simulation,  $ls  are 
estimated  value  of  the  parameters. 


Z(k)  -  ZsimW 
z(k)  -  Zbase 


(11) 


In  Eq.  (11),  Zsim(k)  is  simulation  result. 

The  simulation  precision  is  shown  in  Table  2.  Fig.  7  is  com¬ 
parison  example  between  simulation  result  and  actual  test  data. 
The  validation  results  confirm  that  the  powertrain  simulation 
model  is  reasonably  accurate  and  the  state-space  mathematic 
model  is  suitable  for  math-based  control  system  development. 


3.  Power  control  strategy  development 

Based  on  the  state-space  mathematic  model  shown  in  Eq. 
(9)  and  parameters  estimation  results  shown  in  Table  1 ,  system 
transfer  function  and  pole  points  can  be  figured  out,  as  shown 
below: 


2.4.1.  Sample  data 

According  to  reference  [9],  sample  data  is  very  important  for 
parameters  estimation.  To  make  estimation  result  more  adapted 
to  real  cycle,  in  bench-test  input  command  (namely  input  exci¬ 
tation  signals)  are  designed  elaborately  to  simulate  real  cycle 
condition.  Four  kinds  of  cycle,  including  UDDS,  ECE-EUDC, 
J1015  and  FAST,  are  used  in  bench-test.  Fig.  5  shows  related 
data  in  UDDS  cycle. 

2.4.2.  Parameters  estimation  result 

In  this  work,  the  parameters  are  estimated  based  on  the 
quadratic  least  squares  analysis  by  using  Matlab.  The  estimation 
results  are  shown  in  Table  1 .  To  show  the  estimation  precision 
clearly,  comparison  between  subsystem  model  simulation  result 
and  test  data  is  shown  in  Fig.  6. 


WO)  =  C(sl  -A)~lB 

r  0.126  0  +  0.8819)  0  +  0.1061)  (s  +  0.005547)  -i  T 
5  0  +  0.0534)  0  +  0.73)  0  +  1-25) 

-0.0099181  0  +  0.8819)  (s  +  0.1061)  (s  +  0.005547) 

_  5  0  +  0.0534)  0  +  0.73)  0  +  33.33)  _ 


Table  2 

Model  simulation  precision  (Ee  %)  in  different  cycle 


FAST 

UDDS 

EC_EUDC 

J1015 

mean 

An 

3.91 

2.97 

3.15 

3.31 

3.34 

Alcf 

4.61 

4.32 

3.98 

4.13 

4.26 

4 

6.31 

6.11 

4.73 

5.58 

5.68 

Cbus 

7.74 

6.63 

5.12 

6.13 

6.41 
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time/s 


Fig.  5.  Sample  data  for  parameter  estimation  in  UDDS  cycle. 


si  =0,S2  =  -0.73,  S3  =  -0.0534, 

S4  =  —1.25,  S5  =  —33.33 

In  Eq.  (12)  it  can  be  seen  that  there  is  one  pole  point  equals 
zero.  According  to  Lyapunov’s  stability  theory  [10],  it  can  be 
concluded  that  the  powertrain  is  an  unstable  system.  The  prin¬ 
cipal  task  of  the  control  system  is  to  make  the  system  stable. 
This  is  so-called  power  balance  control  in  practice.  The  power 
control  strategy  must  control  the  battery  SOC  by  influencing  its 
charge/discharge  current  actively,  otherwise  the  SOC  of  the  bat¬ 
tery  may  be  over  upper  limit  or  lower  limit.  This  may  cause  the 
vehicle  anchor  on  the  road. 

Several  academic  research  groups  are  also  working  actively 
on  this  energy  management  technology.  Yuan  Zhu  et  al.  [11] 
solved  the  power  control  problem  of  hybrid  electric  vehicles 
(HEV)  by  fuzzy-based  and  neural  network-based  algorithm. 


This  method  has  good  robustness  and  adaptability,  but  it  needs 
much  more  real-time  operation  and  its  parameters  setting  is 
very  difficulty.  Zaimin  Zhong  et  al.  [12]  proposed  a  battery  cur¬ 
rent  and  SOC  double  close-loop  feedback  algorithm  to  solve 
the  power  balance  control  problem  of  FCHV.  In  this  algorithm, 
the  power  demand  of  electric  motor  was  fed  forward  as  distur¬ 
bance  signal,  and  the  control  parameters  are  set  and  adjusted 
in  vehicle  test.  Though  power  balance  control  is  effective  in 
most  cycle,  this  method  can  be  improved  in  two  respects  at 
least.  Firstly,  Only  the  power  demand  of  electric  motor  was 
fed  back,  others  states  of  system,  such  as  current  of  DC/DC, 
battery  voltage,  and  so  on,  were  not  used  in  control  system. 
Therefore,  it  has  potential  to  optimize  the  control  strategy.  Sec¬ 
ondly,  because  the  algorithm  was  based  on  the  SOC  of  battery 
which  was  estimated  by  integral  of  battery  current,  the  control 
system  was  easy  to  be  affected  by  the  estimated  error  which  was 
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Fig.  6.  Comparison  between  subsystem  model  simulation  and  test  data. 


usually  caused  by  integral  initial  error  and  integral  cumulated 
error. 

Based  on  the  research  results  of  [11,12],  the  paper  pro¬ 
posed  a  power  balance  control  strategy  which  uses  full-states 
feedback  algorithm.  Fig.  8  shows  its  framework.  Where  xjyv. 
control  balance  position,  Uff.  feedforward  command,  Uf^:  feed¬ 
back  command,  x:  states  estimated  by  observer  or  estimator. 
These  variables  have  relationship  shown  in  the  following  equa- 


Fig.  7.  Comparison  example  between  powertrain  model  simulation  and  test  data 
in  UDDS  cycle. 


tion: 

T 

u  —  [  ^fset  fdcfset  Anset  ]  —  uff  +  ufd 

T  T 

Uff  =  [  Pfset_ff  Idcfset-ff  Anset _ff  ]  =  [  0  0  Anset_ff  ] 


i 

0 

\ 

^fsetJd 

0 

fdcfset_fd 

=  K  (xin  —  x)  =  K 

OCVobj 

—  X 

fmset_fd 

0 

V 

0 

/ 

(13) 

where  OCV0bj  •  expect  position  of  the  battery  OCV,  K:  states 
feedback  matrix. 

3.1.  Feedforward  generator 

Feedforward  generator  figures  out  the  feedforward  command 
Uff  and  the  states  balance  position  xju.  In  practice,  Uff  is  figured 
out  according  to  purpose  of  the  driver  basically  and  modified  by 
some  protection  strategy,  such  as  battery  over  discharge/charge 
protection  strategy,  FCE  failure  protection  strategy,  and  so  on. 


Fig.  8.  Framework  of  the  power  balance  control  strategy. 
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time/s 

Fig.  9.  Battery  OCV  estimated  by  the  observer  and  integral. 

One  of  the  state  balance  positions  is  expected  OCV  of  battery. 
Because  SOC  and  OCV  have  corresponding  relationship  in  work 
range  of  the  battery  [7],  the  battery  SOC  can  be  controlled  by 
adjusting  the  state  balance  position. 

3.2.  States  observer 

Because  battery  OCV  and  SOC  have  corresponding  rela¬ 
tionship  in  work  range,  the  battery  SOC  can  be  controlled  by 
feeding  back  battery  OCV.  In  the  state-space  model,  the  state 
U{)  denotes  OCV  of  the  battery.  It  can  be  seen  that  the  state 
U{)  is  very  important  for  control  system.  Unfortunately  the  state 
Uo ,  namely  battery  OCV,  can’t  be  measured  when  powertrain  is 
running.  To  solve  this  problem,  a  Luenberger  state  observer  is 
designed  based  on  the  state- space  model  in  this  paper.  Its  state 
equation  is  showed  in  (14). 

x  =  Ax  +  Bu  +  G(y  —  y)  =  (A  —  GC)x  +  Gy  +  (B  —  D)u 

(14) 

According  to  [10],  no  matter  initial  estimation  error  equals 
zero  or  not,  as  long  as  eigenvalue  of  matrix  A-GC  has  negative 
real  part,  Eq.  (15)  is  always  right. 

lim  [x  —  x]  =  0  (15) 

r— ►oo 

Eq.  (15)  shows  that  if  the  feedback  matrix  of  the  observer 
is  designed  properly,  the  battery  OCV  estimated  by  the  state 
observer  can  tend  to  actual  battery  OCV  by  proper  speed  and 
precision.  Of  course  the  speed  and  precision  is  inconsistent  while 
feedback  matrix  of  the  observer  is  designing.  Because  it  is  inde¬ 
pendent  of  physical  realization  and  the  feedback  matrix  of  the 
observer  can  be  designed  totally  free,  a  most  reasonable  feed¬ 
back  matrix  of  the  observer  always  can  be  designed  finally.  Fig.  9 
is  comparison  between  the  battery  OCV  calculated  by  integral 


Table  3 


Signification  of  the  state  feedback 


Feedback 

Signification 

K\  1  (/<dcf  —  4i)  ^21  (7dcf  —  4n) 

Compensating  delay  of  7dcf  caused 
by  FCE  slow  responding  feature 

/m 

Making  FCE  power  set  following 
load  power  demand 

Ki3Uo,K23Uo 

Control  SOC  of  the  battery 

KuUuKuUi 

Compensating  the  slow  part  of  the 
battery  voltage 

K\5  u2,  k25  u2 

Compensating  the  slow  part  of  the 
battery  voltage 

and  that  of  estimated  by  the  observer  in  a  J1 01 5  cycle  vehicle 
test.  According  to  [7],  during  short  time  range,  the  battery  OCV 
calculated  by  integral  equals  to  the  actual  value  approximately. 

It  can  be  seen  from  Fig.  9  that  the  OCV  estimated  by  observer 
tend  to  actual  battery  OCV  by  moderate  speed  and  precision. 
It  shows  that  the  feedback  matrix  of  the  observer  designed  is 
reasonable. 

3.3.  States  feedback  controller 

The  key  of  states  feedback  controller  is  the  states  feedback 
matrix,  as  shown  in  Eq.  (13).  The  states  feedback  matrix  is  a  three 
rows  and  five  columns  matrix,  and  each  element  of  the  matrix 
has  special  signification  which  is  shown  in  Table  3.  Where,  Kx j 
is  element  of  the  states  feedback  matrix. 

The  states  feedback  matrix  can  be  designed  without  any  con¬ 
straint  in  control  theory,  but  in  practice,  it  needs  to  be  constrained 
by  actual  conditions  which  are  showed  in  Eq.  (16). 

*11  -  *12  =  *21  -  *22  =  1 
*13  =  *23  <  a 

*14  =  *24,  *15  =  *25  <  b  (16) 

*11  <  *21  <  C 

*31  =  *32  =  *33  =  *34  =  *35  =  0 

where  a,b,c:  experimental  upper  limit  determined  by  actual  con¬ 
dition.  The  fifth  equation  of  the  Eq.  (16)  means  that  the  electric 
motor  torque  is  set  according  to  feedforward,  as  shown  in  Eq. 
(13).  With  the  constraints  defined  in  Eq.  (16),  each  element  of 
the  matrix  is  selected  respectively.  Proper  K\\  and  K21  can  com- 
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Fig.  10.  Step-responding  curve  of  the  powertrain  in  difference  A3  3  and  K23. 
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Fig.  11.  Time-trace  of  vehicle  speed  (UDDS  cycle). 


Fig.  13.  Time-trace  of  battery  SOC  (UDDS  cycle). 


pensate  delay  of  7dcf  caused  by  FCE  slow-responding  feature 
and  cut  FCE-responding  time.  But  over  large  K\\  and  K2\  may 
cause  FCE  power  setting  vibrating  in  some  condition.  In  prac¬ 
tice  K\\  and  A^i  are  selected  according  to  the  experiment.  K\^ 
and  7^23  are  very  important  for  state  feedback  controller  because 
they  determine  the  stability  and  stable  level  of  the  control  sys¬ 
tem.  By  analyzing  the  step  responding  of  the  powertrain,  K\^  and 
K22  can  be  selected  properly.  Selected  system  step-responding 
curve  are  presented  in  Fig.  10.  It  can  be  seen  that  the  bigger  K\2 
and  K22  are,  the  faster  converging  behavior  of  Uq  is.  But  over 
big  K\2  and  K22  will  cause  battery  voltage  vibrating.  As  shown 
in  Fig.  10,  when  K\2  and  K22  equal  10,  converging  speed  of 
U{)  is  moderate  and  battery  voltage  converges  to  zero  simplex. 
Like  K\ 3  and  K22,  K 14, ^15,^24,^25  can  be  selected  by  the  same 
way. 

4.  Validation  the  power  control  strategy  by  field  test 

Based  on  the  power  balance  control  strategy  designed  above, 
the  powertrain  control  strategy  is  developed  consequently  and 
applied  in  START3  prototype  vehicle.  To  validate  the  strategy, 
simulation  and  field  test  are  carried  out  using  different  cycle, 
including  J1015,  UDDS,  ECE-EUDC.  The  result  of  the  simu¬ 
lation  and  field  test  shows  that  the  powertrain  control  system  is 
stable  and  battery  SOC  fluctuates  between  ideal  upper  limit  and 
lower  limit.  Selected  simulation  and  test  data  of  UDDS  cycle 
are  presented  in  Figs.  11-13. 

Example  time-trace  of  the  subsystems  current  in  simulation 
and  field  test  are  shown  in  Fig.  12.  It  can  be  seen  that  the  fuel 
cell  provides  the  bulk  of  the  power  load  while  the  battery  helps 
with  the  transient  part.  Fig.  13  shows  that  battery  SOC  in  UDDS 
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Fig.  12.  Time-trace  of  sub-systems  current  (part  of  the  UDDS  cycle). 


cycle  fluctuates  between  48%  and  52%.  The  object  of  powertrain 
control  system  is  achieved  excellently. 

5.  Conclusions 

A  dynamic  state- space  mathematic  model  of  FCHV  pow¬ 
ertrain  is  developed  based  on  the  lab-testing  and  field  vehicle 
testing  result  from  START3  prototype  FCHV.  A  vehicle  simu¬ 
lation  model  is  constructed  by  following  the  configuration  used 
in  START3.  Model  validation  results  confirm  that  the  simula¬ 
tion  model  is  reasonably  accurate  and  the  mathematic  model  is 
suitable  for  model-based  control  strategy  development. 

Based  on  the  state- space  model,  a  power  balance  control 
strategy  is  developed,  which  uses  full-states  feedback  algorithm. 
A  Luenberger  state  observer  is  designed  to  estimate  the  OCV 
of  the  battery,  which  makes  state-feedback  feasible.  The  state 
feedback  controller  is  then  designed.  Results  of  simulation  and 
vehicle  field  test  show  that  performance  of  the  control  system  is 
excellent. 

The  process  described  in  this  paper  forms  the  model-based 
FCHV  powertrain  control  system  development  method,  which 
is  applied  in  development  of  START3  FCHV  powertrain  con¬ 
trol  system  successfully,  and  which  is  the  general  guidance  for 
powertrain  control  strategy  development  for  FCHV. 

Theoretically  full-states-feedback  method  is  advanced, 
which  enhances  degree  of  freedom  in  control  system  design,  but 
robustness  of  the  control  system  is  ordinary,  and  in  this  paper, 
index  of  control  system  performance  is  mainly  discussed  in  fre¬ 
quency  domain.  These  problems  will  be  discussed  in  subsequent 
research  project. 
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